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The entanglement properties of two-electron atomic systems have been the subject of considerable 
research activity in recent years. These studies are still somewhat fragmentary, focusing on numerical 
computations on particular states of systems such as Helium, or on analytical studies of model- 
systems such as the Moshinsky atom. Some general trends are beginning to emerge from these 
" studies: the amount of entanglement tends to increase with energy and, in the case of excited 

states, entanglement does not necessarily tend to zero in the limit of vanishing interaction between 
. the two constituting particles. A physical explanation of these properties, shared by the different 

■ two-electrons models investigated so far, is still lacking. As a first step towards this goal we perform 

here, via a perturbative approach, an analysis of entanglement in two-electrons models that sheds 
new light on the physical origin of the aforementioned features and on their universal character. 
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I. INTRODUCTION 

^: 

Entanglement constitutes one of the most fundamental phenomena in Quantum Mechanics [H-Q • Entangled states 
of composite quantum systems exhibit non-classical correlations that give rise to a rich variety of physical phenomena 
■ of both fundamental and technological significance. Quantum entanglement can be considered in two different and 
^ | complementary ways. On the one hand, entanglement can be viewed as a resource. The controlled manipulation 
O" 1 of entangled states is at the basis of several quantum information technologies. On the other hand, entanglement 
can be regarded as a fundamental ingredient for the physical characterization of natural quantum systems such as, 
I . for instance, atoms and molecules (for a comprehensive and up to date review on this subject see [fj). These two 
1 points of view are closely related to each other, although the latter is somehow less developed than the former. 
ON . Concerning the second of the abovementioned approaches, several researchers have investigated in recent years the 
ON ■ phenomenon of entanglement in atomic physics 6H17| . This line of enquiry is contained within the more general 
program of applying tools and concepts from information theory to the analysis of atomic and molecular systems [l8l - 



27J . Most of the studies on entanglement in two-electron systems focused on the properties of the concomitant ground 
ON ' states. However, the entanglement features exhibited by excited states of two-electron atomic systems have also been 
explored [HI, [l4[ . In this regard, the most detailed results have been obtained from analytical investigations of the 
entanglement properties of exactly soluble models, specially the Moshinsky one [l3|. The behaviour of these soluble 
models is consistent with some (partial) results yielded by numerical explorations of entanglement in Helium based on 
high quality, state-of-the-art wave functions. Some general trends begin to emerge from these investigations. First, 
and not surprisingly, entanglement is found to increase with the strength of the inter-particle interaction. Second, 
entanglement also tends to increase with energy. Finally, the entanglement of excited states does not necessarily 
vanish in the limit of zero interaction. The last two properties are, perhaps, less intuitively clear than the first one. 
In fact, in a recent comprehensive review article on entanglement in atomic and molecular systems by Tichy et al. @ 
it is said that ". . . the limit of vanishing interaction strengths does not necessarily yield a non-entangled state . . . it 
remains open whether this discontinuity effect has to be considered an artefact of the entanglement measures that are 
used, or whether a physical explanation will be provided in future." . 

It is remarkable that the various two-electron models where entanglement has been studied so far share the basic 
qualitative features mentioned above. This suggests that these features may constitutegeneric properties of these kind 
of two-fermion models. In this regard, we share the opinion expressed by Tichy et al. [6[: "Due to the non-integrability 
of any non-hydrogen-like atom, theoretical studies of multielectron systems have, so far, mainly focused on exactly 
solvable model atoms. While such models differ strongly from real multielectron atoms as concerns the interelectronic 
interaction and the definition of the confining potential, they allow insight in some qualitative features." 

The aim of the present work is to clarify the origin of the aforementioned properties. To this end we are going to 
consider a perturbative approach to this problem, regarding the term in the Hamiltonian describing the interaction 
between the two electrons as a small perturbation. We shall show that the eigenvalue degeneracy of the unperturbed 
Hamiltonian (describing independent particles) plays a crucial role in explaining the entanglement features of the 
"perturbed" system. 
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II. QUANTUM ENTANGLEMENT IN SYSTEMS OF TWO IDENTICAL FERMIONS 

Correlations between two identical fermions that are only due to the antisymmetric nature of the two-particle state 
do not contribute to the state's entanglement [28 36]. The entanglement of the two-fermion state is given by the 
quantum correlations existing on top of these minimum ones. For example, a two-fcrmions state of Slater rank one 
(that is, a state whose wave function can be expressed, in terms of an appropriate singlc-particc orthonormal basis, as 
one single Slater determinant) must be regarded as non-entangled. There are deep, fundamental physical reasons for 
this. On the one hand, the correlations exhibited by such states are not useful as a resource to perform non-classical 
information transmission or information processing tasks [28] . On the other hand, the non-entangled character of 
these states is consistent with the possibility of associating complete sets of properties to both parts of the composite 
system (see [29l - l3ll | for a detailed analysis of various aspects of this approach). 

Two useful quantitative measures for the amount of entanglement of a pure state of a system of two identical 
fermions are expressed in terms of (see (37j and references therein) the linear entropy, 

e L {\i>)) = l-2Tr\p% (1) 

and the von Neumann entropy 

s vN (\ip)) = -Tr\p r lnp r }-ln2, (2) 

of the single-particle reduced density matrix p r . Notice that according to the entanglement measures given by Eqs. 
(TI]) and ^ a pure state that can be represented by a single Slater determinant has no entanglement (that is, it is 
separable). 

The fermionic entanglement measures ([1]) and J2| are closely related to the Schmidt decomposition of pure states 
of systems constituted by two identical fermions |32j . For any pure state \ip) of two identical fermions it is posible to 
find an orthonormal basis i = 0, 1, . . .} of the single-particle Hilbert space such that the state l^) can be written 



E 



(|2<>|2i + l)-|2i+l)|2i>), (3) 



where the Schmidt coefficients Aj verify < Aj < 1 and J^. Aj = 1 (in the case of systems with a single-particle Hilbert 
space of finite dimension TV, we assume that N is even and that the sums on the index i run from i = to i = N/2). 
Then one has that the entanglement measures ([I]) and ([2]) can be expressed in terms of the Schmidt coefficients of 
the state \tp) respectively as (H, HtJ, 

Sl(\4>)) = l-J2 X l ( 4 ) 

i 

and 

£vn(\tP)) = ~y^A»lnAi. (5) 

i 

In the particular case of systems of two fermions with a single-particle Hilbert space of dimension four, the quantity 
2el reduces to the entanglement measure (usually referred to as squared concurrence) studied in [28| (see also [3o|). 
The entanglement measure given by equations ((TJ) and (jU has been recently applied to the analysis of various physical 



systems or processes, including electron-electron scattering processes [34[, the study of entanglement-related aspects 
of quantum brachistochrone evolutions [35j , and the entanglement properties of two-electron atomic models [HI . As a 
final remark on entanglement in fermionic system we mention that in the present work we deal with the fermionic case 
of the concept of entanglement between particles. This is not the only possible conception of entanglement in systems 
of identical particles. In particular, there is an approach to the study of entanglement in systems of indistinguishable 
particles which focuses on the entanglement between different modes (see, for instance, [6j and references therein). 



III. PERTURB ATIVE APPROACH 



Let us consider a system of two identical fermions ("electrons") governed by a Hamiltonian of the form H = 
Hq + AH' , where the unperturbed Hamiltonian Hq corresponds to two independent (non-interacting) particles, AH' 
describes the interaction between the electrons, and A is a small parameter. When this system is treated perturbatively, 
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the pcrturbative corrections to the eigenenergies correspond to some "fine structure" sitting on top of the main pattern 
due to the spectrum of H . It is plain that within this scenario the leading, zeroth-order contribution to the energy 
spectrum is independent of the detailed structure of the perturbation H' . As we shall presently see, the situation 
is completely different when, instead the energy, we calculate the entanglement of the system's eigenstates. When 
the unperturbed energy eigenvalues are degenerate the leading (zeroth-order) contribution to the eigenfunction's 
entanglement does depend, in general, on the details of the perturbation. 

Let us consider an m-fold degenerate energy level of Hq, with an associated set of m orthonormal eigenstates 
\tpj), j = 1, . . . m. Since Hq describes two non-interacting particles, the m eigenstates \tpj) can always be chosen 

( 1 2) 

to be Slater determinants written in terms of a family of orthonormal single-particle states \<j)j ' ), so that = 

( 1 / a/2 ) ( | ^* ) I ) — l^j 2 ^)!^ 1 ^))- All the members of the subspace H s spanned by the states \tpj) are eigenstates 
of Hq corresponding to the same eigenenergy. That is, energywise they are all equivalent. However, the different 
members of this subspace have, in general, different amounts of entanglement. Typically, the interaction H' will lift 
the degeneracy of the degenerate energy level. If we solve the eigenvalue problem corresponding to the (perturbed) 
Hamiltonian H and take the limit A — > 0, the perturbation H' will "choose" one particular basis {|V>j.)a->-o} among 
the infinite possible basis of W s . The states constituting this special basis will in general be entangled. These states 
are of the form \i/j' k )x^-o = Ej=i c kj\ipj), an d ar e determined (according to standard perturbation theory [38|) by the 
eigenvectors of the m x m matrix H with elements given by Hij = (ipi\H'\ipj) . It is then clear that in the limit A — > 
the eigenstates of H will in general be entangled. 

( 1 2) 

Let m be the number of different single-particle states within the family {|<^ ' ; ), 1, . . . , to}. It is a quite typical 
behavior that to tends to increase with the degree of degeneracy to of the energy levels of Hq which, in turn, tends 
to increase with energy (that is, it tends to increase as one considers higher excited states). This explains (at least 
in part) why the range of entanglement-values available to the eigenstates {|?/>JJ.)a->o} tends to increase with energy. 
Indeed, the maximum amount of entanglement (as measured by @) that can be achieved by a linear combination of 

(1 2) 

Slater determinants constructed from the single-particle states ' ')} is 

e = lnfi, (6) 

where fl is the integer part of to/2. Expression (jSJ) provides an upper bound for the entanglement of the states 
M)a^o}. 

In the present work we are going to focus on the entanglement properties exhibited by the eigenstates {|^.)a-s-o} 
and on the entanglement upper bound ([6]). In this regard, our pcrturbative approach is unusual, since we are focusing 
on "zeroth-order properties". However, it must be stressed that the amounts of entanglement of the states {IV^a^o} 
are in general finite quantities that do not vanish when A — > and, consequently, constitute dominant aspects of the 
entanglement-related features characterizing the system. 



IV. TWO INTERACTING SPIN-i FERMIONS IN AN EXTERNAL CONFINING POTENTIAL 

We apply now the previous considerations to a system consisting of two interacting spin-^ fermions in an external 
confining potential U(x). The interaction between the particles is described by the potential function V(x% — x%), 
with V an even function. The Hamiltonian of this system is then, 

H = ~§ J? " \ih + u{xi) + U{X2) + xv(xi ~ X2) ' (7) 

where x\ and Xi are the coordinates of the two particles. We use atomic units (m = 1, h = 1). A relevant instance 
of this system corresponds to the case of harmonic confinement, U(x) = h^x\ where w is the natural frequency of 
the external harmonic field. This case includes the Mishonsky atom 0, UH, |39| , where the interaction between the 
particles is also harmonic, \V(x\ — x^) = |Aw 2 (a;i — X2) 2 , with Alu 2 > being the square of the natural frequency of 
the interaction harmonic field. The Moshinsky atom is an exactly soluble system whose entanglement properties have 
been studied in detail. The examples considered here indicate that some important entanglement-related features of 
the Moshinsky model are also encountered in more general systems. 

We now apply the formalism of perturbation theory to a system described by the Hamiltonian (|7J| with harmonic 
confinement and a generic interaction V between the particles. The unperturbed Hamiltonian is then, 

1 d 2 1 d 2 1 2 2 1 2 2 
Fo = -2^-2^ + 2^ + 2^' (8) 
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and the perturbation, 

XH' = XV(x 1 -x 2 ). (9) 

When A = 0, the model consists of two independent harmonic oscillators with the same natural frequency. Let \n) (n = 
0, 1, 2, ...) be the eigenstatcs of each of these oscillators. Then, the kets \n, ±) constitute a single-particle orthonormal 
basis (the signs ± correspond, in standard notation, to the spin state of the spin-^ particle). The eigenstates of Hq 
are characterized by two quantum numbers n\ and n 2 corresponding to the alluded pair of independent oscillators. 
The corresponding eigenenergies depend only on the value of the sum n\ + n 2 and are m-fold degenerate with 
m = 2{n\ + n 2 ) + 1 (m = 2(ni + n 2 ) + 2) if ri\ + n 2 is even (odd). Assuming that n\ + n 2 is odd, with ri\ = n 2 — 1, 
and taking spin into consideration, we can choose the following set of m antisymmetric eigenstates (all with the same 
energy), 

IV>i>= 75 (|»i,+)|na,+)-|n a ,+)|ni J +)) 
|V'2>= (\n 1 ,+)\n 2 ,-}-\n 2 ,-}\n 1 ,+)) 
H's)= -= (|ni,-)|n2,+)-|n2,+)|ni,-» 
1^4)= 7= (|ni,-)|n 2 ,-) - |n 2 ,-)K,-)) 

i (10) 

\1> m -s)= 75 (|0,+)|n 2 + n 1 ,+)-|n 2 + ni,+)|0,+)) 

|Vm-2>= -= (|0,+)|n 2 + ni,-)-|n 2 + m,-)|0,+)) 

|Vm-i)= 7f (|0 ) -)|n 2 + m,+)-|n 2 + ni,+)|0,-)) 

\i>m)= -= (|0 ! -)|n 2 + n 1 ,-)-|n 2 + m,-)|0,-)), 

which are represented by single Slater determinants and consequently have zero entanglement. A similar set of 
separable eigenstates of Hq can be chosen when n\ + n 2 is even. 

We consider now a harmonically confined two-fermion system with an interaction potential given by a repulsive 
Dirac delta function, 

XV = ^X5( Xl - x 2 ), (11) 

For the first excited energy level of Hq (ni + n 2 = 1) which is four-fold degenerate we then have, 

_ /0 0\ 

Vo 00/ 

and the corresponding eigenvectors can be written as 

= IV> 4 ) 
1^3) = IV-i) 

|^) = -L(|^> + |^ 3 )). (13) 

In the limit of vanishing interaction, A — > 0, the eigenstates corresponding to the first two excited energy levels of the 
full Hamiltonian tend to the states (fT3"|) . which have the following amounts of entanglement, 

eM))=s L (m) = (14) 
£vN(m)=e vN (\ip' 3 ))=0 (15) 

e L {W 1 ))=e L (W i )) = \ (16) 
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vN (\tp[)) = e„Jv(M» = ln2 



(17) 



It can be verified after some algebra that the eigenvectors (and the associated amounts of entanglement) obtained in 
this case coincide with those corresponding (in the A — > limit of the two first excited energy levels) to an harmonic 
interaction. That is, they are the same as those associated with the Moshinsky model. We thus see that the harmonic 
and the Dirac delta interactions lead, for particles confined by an external harmonic well, to the same entanglement 
behaviour of the first excited states in the limit of weak interaction. 

We now consider the case of generic external (one dimensional) potential U(x) and interaction V{x\ — X2) (with V 

an even function) so that Hq = Ym=i fjj? + ^O^i)) an d A-ff' = W(xi — X2), which leads to 



H 



fa 
b c 
c b 

\0 







(18) 



with 



a = -((0|<l|i/'|0)|l)-<0|<l|i/'|l)|0)-<l|(0|i/'|0)|l) 
b = i((0|(l|H'|0)|l) + (l|(0|ff'|l)|0)) 
c = -i((0|(l|i/'|l)|0) + (l|(0|iJ'|0)|l)) 



(l|(0|ff'|l)|0)) 



(19) 



where |0) and |1) are the ground and first excited eigenstates corresponding to the external confining potential U(x) 
The eigenvectors of H are, 



1 

71 

1^4) 

IV-i) 

T2 m 



1^3)). 



(20) 



The values of entanglement exhibited by these states are Sl = £ v n = associated to \ip' 2 ) and 1^3) and sl = \, 
e V N = In 2 associated to |^) and \ip'^). This result generalizes the previous one, as we have solved the problem for 
generic interactions and external confining potentials. Here it is possible to obtain general results for an arbitrary 
confining potential U (x) for the case where one has (in the limit of vanishing interaction) one particle in the ground 
state and one particle in the first excited state of U(x), because the degeneracy of the concomitant energy level (of 
the two-partile system) can be determined directly without knowing the detailed energy spectrum associated with 
U(x). On the other hand, the properties exhibited by states of higher excitation in the limit of vanishing interaction 
do depend (via the degeneracy appearing in this limit case) on the detailed eigenenergies of the confining potential 
U(x). Consequently, the analysis of the limit of vanishing interaction can be performed only in a case-by-case way. In 
the next section we are going to consider higher excited states in the case of a generic interaction potential V{x\ — X2) 
and a harmonic confining potetial. 



V. ENTANGLEMENT UPPER BOUND FOR EXCITED STATES IN THE LIMIT OF WEAK 

INTERACTION 

We consider now two spin-^ particles (in one dimension) confined by an external harmonic potential and having a 
generic interaction W(xx — x 2 )- We shall calculate general upper bounds for the entanglement of the eigenstates of 
this system in the limit A — > 0. These bounds, expressed in terms of the quantum numbers ni and n 2 characterizing 
the eigcnf unctions of Hq, are, 

e L {\n in2 )) < 1 2 (21) 
n\ + 7i 2 + 1 



£„jv(|nin2)) < ln(ni + n 2 + 1) 



(22) 
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Equation (|22[) constitutes a particular instance, corresponding to a harmonic confining potential, of the general upper 
bound ©. In Fig. Q]we plot the entanglement bounds against ni+n 2 - These curves represent the maximum possible 
entanglement compatible with those quantum numbers (the bounds do not depend on the interaction and are, in this 
sense, universal). 




FIG. 1: Entanglement bound as a function of the sum of the quantum numbers 711+712. All depicted quantities are dimensionless. 



In the particular case rii+712 = 2, besides the above upper bound, we can calculate the exact amount of entanglement 
(in the limit A — > 0) for an arbitrary interaction potential V{x\ — x 2 )- This case corresponds to a five- fold degenerate 
energy level of Hq shared by a set of eigenvectors of the form (|10p . with m = 5. The generic matrix H is of the form 



with 



(a \ 

b c d ' 

H = c b -d 

a 

\0 d -d e / 



a =-((01(21^10)12) - (01(21^12)10) - (2|(0|H'|0)|2) + <2|<0|iJ'|2)|0)) 



(23) 



b =-((01(21^10)12) + <2|<0|tf'|2)|0)) 
c = -l((0|(2| J ff'|2)|0) + (21(01^10)12)) 



(24) 



d=-((0K2|H'|l)|l) + (2K0|H'|l)|l)) 
e=(l|(l|ff'|l)|l). 

The corresponding normalized eingevectors can be expressed as follows 

Wi) = IV>1> 



W2) 



1^4) 

^(i^) + i^)) 



1 



1 



1 



(-ri|V>2>+n|V>3> + |V>5» 



(-r 2 |V>2) +r 2 \i/> 3 ) + ^5)), 



(25) 
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where r\ and r 2 are functions of the matrix elements b, c, d, e given by the eexprcssions, 

C.+cR d-cR 
ri = C^TdR' r2 = C^dR' (26) 

with 

d = -be + c 2 + 2d 2 + ce, C 2 = d{b + 3c-e), (27) 

and 



R = ^b 2 + c 2 + 8d 2 + 2ce + e 2 - 2b(c + e) (28) 

Now we calculate the amounts of entanglement of the states l^) and \ip' 5 ) that depend (in the same way) on T\ and 
ri respectively. The values adopted by these two constants depend on the form of the interaction V(x\ — x-i). The 
general expression for the amount of entanglement of the states {ip'^} and are, 

£L{n) = l -j0[Th 1=1,2 (29) 



4r 2 

e vN {n) = ln(2r? + 1) - lnr, i = 1, 2 (30) 

We plot these expressions in Fig. [2j Note that the particular value r; = corresponds to the harmonic interaction 
in the Moshinky model. 




FIG. 2: el (left panel) and e v n (right panel) as a function of for the case n\ +ri2 = 2 and an arbitrary interaction V(xi — £2). 
The stars in both plots correspond to the entanglement amount for the Moshinsky atom (harmonic interaction). All depicted 
quantities are dimensionless. 



VI. CONCLUSIONS 



By recourse to a perturbative approach we studied the entanglement-related properties of a system consisting of 
two interacting spin-^ fermions ("electrons") confined by an external potential. Our present results clarify some 
aspects of the basic entanglement features exhibited by particular two-electrons models studied previously, and 
shed some light upon the fact that these systems share important qualitative entanglement properties that are also 
observed in more general cases. Our analysis highlights the important role played by the degeneracy of the energy 
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levels of the "unperturbed" (interaction- free) Hamiltonian Hq. The non- vanishing entanglement exhibited by the 
interacting particles in the limit of vanishing interaction is due to the particular eigenbasis of H "chosen" by the 
interaction. This amount of entanglement tends to increase with the alluded degeneracy which, in turn, tends 
to increase with energy. This sheds light on the physical reasons behind the fact (observed in all cases studied so 
far) that the amount of entanglement exhibited by the eigenstates of two-electrons systems tends to increase with 
energy. These basic trends do not depend on the particular entanglement measure employed, as has been shown 
in the present work, where entanglement measures based upon the linear and the von Neumann entropies were 
considered. In connection with this point it is worth to mention a relevant question raised in a recent review article 
on entanglement in atoms and molecules @: does the existence, for some excited states, of a finite amount of 
entanglement in the limit of vanishing interaction depend upon the particular entanglement measure employed? As 
already mentioned, the answer to this question is that the alluded feature constitutes an intrinsic property of the 
systems under consideration, which does not depend on the entanglement measure used. 

As particular illustrations of the above considerations we have 

• Computed for general confining and interaction potentials, in the limit A — > 0, the entanglement measures 
based upon the linear entropy and the von Numann entropy of the excited states associated with the four-fold 
degenerate unperturbed first excited state. 

• Obtained for a harmonic confining potential and a generic interaction potential V{x\ — X2), the entanglement 
of the eigenstates corresponding in the above limit to the second excited unperturbed energy level (given by 
n\ + ri2 = 2). 

• Determined, for a harmonic confining potential and an arbitrary interaction, upper bounds on the amounts of 
entanglement exhibited by the system's eigenstates in the limit of vanishing interaction. These upper bounds 
are expressed in terms of the quantum numbers n\ and «2- 

The results advanced here corresponding to the entanglement in the limit of vanishing interaction are exact, and 
our procedure can in principle be applied to any excited state of systems of the kind considered in this work. In this 
sense, the perturbative method used here is not (in intsclf) the fundamental protagonist of our present considerations. 
The perturbative approach was used only as a tool to determine (exactly) the entanglement features of the system's 
eigenstates in the limit of vanishing interaction, in order to get some insight on the qualitative entanglement features 
of two-elctrons models. It is worth stressing that the entanglement exhibited by these systems in the limit A — > 
constitutes a basic, dominant aspect of their entanglement-related characteristics. It is to be expected, on general 
physical grounds, that the entanglement-degeneracy relationship uncovered here constitutes a typical feature of atomic 
models. This provides a first step towards explaining the fact that the general entanglement features exhibited by 
soluble models such as the Moshinsky one are also observed in other two-electrons systems. 

The ideas discussed here may also be useful for the analysis of entanglement-related aspects of other scenarios 
involving interacting fermions, such as those appearing in molecular or solid state physics. Just to mention one 
example, a situation similar to the one observed in the atomic models considered here also occurs when studying the 



behaviour of electronic entanglement in the dissociation process of diatomic molecules 27j. One observes that, for 
instance, in the limit of large values of the reaction coordinate (corresponding to vanishing interaction between the 
electrons in the system) describing the dissociation of H2, the electronic entanglement does not tend to zero f27j . 
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